Bond breaking and bond formation: how electron correlation is captured in 
many-body perturbation theory and density-functional theory 
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For the paradigmatic case of Jb-dissociation we compare state-of-the-art many-body perturba- 
tion theory (MBPT) in the GW approximation and density-functional theory (DFT) in the exact- 
exchange plus random-phase approximation for the correlation energy (EX+cRPA). For an unbiased 
comparison and to prevent spurious starting point effects both approaches are iterated to full self- 
consistency (i.e. sc-RPA and sc-GW). The exchange-correlation diagrams in both approaches are 
topologically identical, but in sc-RPA they are evaluated with non-interacting and in sc-GW with 
interacting Green functions. This has a profound consequence for the dissociation region, where sc- 
RPA is superior to sc-GW\ We argue that for a given diagrammatic expansion, the DFT framework 
outperforms the many-body framework when it comes to bond-breaking. We attribute this to the 
difference in the correlation energy rather than the treatment of the kinetic energy. 
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First-principles electronic-structure calculations have 
become indispensable in many fields of science, be- 
cause they yield atomistic insight and are complemen- 
tary to purely experimental studies. Since the full 
many-body problem of interacting electrons and nuclei 
is intractable for all but the simplest systems, different 
strategies for approximate approaches have been devel- 
oped over the years. The most prominent are density- 
functional theory (DFT) 0-0], many-body perturbation 
theory (MBPT) coupled-cluster theory Q and 

quantum Monte Carlo methods ||. Each approach has 
its strengths and weaknesses in terms of accuracy, appli- 
cability, and computational efficiency and no consensus 
has been reached regarding the optimal approach for cur- 
rent and future challenges in electronic-structure theory. 
In this work we address the difference between DFT and 
MBPT for the total energy and ask the questions: Given 
a fixed set of diagrams for the electron-electron interac- 
tion, will the DFT and the MBPT framework give the 
same result? And if not, which one is better? To answer 
these questions we consider the paradigmatic case of H 2 
dissociation. 

In the past, DFT and MBPT have been compared 
directly in the exchange-only case 11|. In MBPT this 
corresponds to the Hartrcc-Fock approach, whereas in 
DFT a multiplicative Kohn-Sham (KS) potential is con- 
structed by means of the optimized effective potential ap- 
proach (OEP) [12 ]. However, this comparison is mostly 
of academic interest, because electron correlation is not 
included. As we will demonstrate in this Letter, the 
comparison between DFT and MBPT can be extended 
to encompass correlation using exact-exchange plus cor- 
relation in the random-phase approximation to DFT 
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Figure 1. $ functional for RPA and GW correlation ener- 
gies (Eq. [SJ. The arrowed lines correspond to the interacting 
Green function G in GW , and the KS Green function G s 
in RPA. Dashed lines denote the bare Coulomb interaction, 
and the minus sign of the prefactor comes from the rules for 
evaluating Feynman diagrams [9. Ilfj]. 



(EX+cRPA), [TKll] referred to as RPA in the following, 
and the GW approach to MBPT @,[l7j]. The exchange- 
correlation diagrams in both approaches are topologically 
identical (see Fig. Q}, but in RPA they are evaluated 
with a non-interacting KS and in GW with an interacting 
Green function. To illustrate the impact of these differ- 
ences we consider the bond-breaking/formation regimes 
in the binding curves of H2. To avoid starting point 
effects both approaches are iterated to self-consistency, 
which we denote as sc-RPA and sc-GW. 



Let us start with the ground-state total-energy expres- 
sion for an interacting electron system obtained with the 
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Figure 2. (Color online) Total energy (eV) of the H2 molecule as a function of bond length (A). Different flavors of GW 
and RPA are shown compared to PBE, rPT2 and accurate full configuration interaction calculations taken from Ref. fl8l ]. 
Hartree-Fock (HF) and exact-exchange OEP (OEPx) are identical for H2 and are included for comparison. All calculations 
were performed using a Gaussian cc-pVQZ [19( basis set. 



adiabatic-connection (AC) technique (see e.g. Ref. [la]): 



E =E - ^J d\ J drdr'v(r,r')x 

[ — XA(r,r';iw) + n(r)5(r-r') 
Jo n 
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Here v(r,r') is the Coulomb interaction, Tr [AB] denotes 
/ drdr'A(r,r')B(r',r) and E Q = T s + E H + E cxt , T s being 
the kinetic energy of the KS independent-particle system, 
En the Hartree and E ext the external energy. Along the 
AC path (i.e. at each value of A), the electron density 
n(r) is assumed to be fixed at its physical value, and 
Xai G\, and are the polarizability, the single-particle 
Green function, and the self-energy, respectively. 

The RPA for the total energy can be most conve- 
niently introduced in Eq. [T] through the approximation 
X\ = Xs(l - Avxs)^ 1 , where \s = X\=o = -iG s G s , with 
G s being the non-interacting KS Green function also of- 
ten denoted as Go- Within this approximation, the inte- 
grand in Eq. [T] is assumed to depend on A only through 
the scaled Coulomb interaction Aw. Alternatively the 
RPA total energy can also be obtained through Eq. [2] 
by introdu cing the GW approximation for the proper 
self- energy 16. 2fj|: 



du' 
27 



G A ( W + W ')W A ( W ')e 1 "'", (3) 



where Wa[Ga] = Aw(l + iXvGxGx) -1 and rj is a posi- 
tive infinitesimal. The RPA total energy is retrieved by 
omitting the A-dependence of G\, i.e. replacing G\ by 
the KS non-interacting Green function G s = Ga=o and 
W\ by W s = W\[G S }. Either way, the A integration in 
Eqs. Q] or [2] can now be carried out, yielding the sum of 
the exact-exchange energy E x and the RPA correlation 



energy Ef PA , where: 

Ef PA = r^Tr[ln(l 
Jo 2lT 



Xs(iu)v) + Xs{iu)v] . (4) 



Combining Eqs. Q] and [4] allows us to express the RPA 
total energy functional as: 



E [G s ] =T S + E cxt + En + E K 



E, 



RPA 



(5) 



We now come to the differences in the evaluation of the 
total energy in the context of KS-DFT and MBPT. In 
MBPT, the Green function G\ represents an interacting 
electron system for A ^ 0, and has to satisfy the Dyson 
equation: 



g: 



Ea[G> 



+ "ext + (1 - A)l)H + V xc (6) 



with w A xt being the external potential of the A-dependent 
system (chosen to keep the density fixed), and v xc the 
exchange-correlation potential of the KS non-interacting 
particle reference system. Making use of Eqs. [3] and [6j 
the A-integration in Eq. [2] can be carried out and one 
arrives at the following expression for the total energy 

E = -Eg[G\ + $[G] - — / dujx 

Tr [(G^M + v xc )G(iuj) - 1 + ln(G- l (io;))] . (7) 

Details for the derivation of Eq. [7] can be found in the 
supplemental material [2l[. In Eq.fTl G = G\ = i, and the 
functional $ [G] is defined as 0, lioj ] 
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where S( n ' is the sum of all self-energy diagrams that 
contain n explicit Coulomb interaction lines. We note 
that since S = S&/5G , an approximation for $ directly 
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Table I. Total energies (in eV) of H2 in the equilibrium geometry and deviation (A) from the exact reference. 
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translates into a corresponding approximation for S. The 
diagrammatic representation of $ in the GW approxima- 
tion is illustrated in Fig. Q] 

In the KS framework, the sc-RPA total energy is ob- 
tained by requiring G s in Eq.[S]to satisfy the Dyson equa- 
tion G s (iui) = (iu> + V 2 /2 — v CKt — vh 



,RPA\ 



where 



,RPA 



is determined by the optimized effective potential 
equation (als o known as the linearized Sham-Schliitcr 
equation) jH,^: 



G S (E GW [G S ] 



,RPA 



)G S = 0. 



(9) 



Alternatively, the sc-RPA energy can be obtained by min- 
imizing £ ,RPA in Eq.[5]with respect to the non-interacting 
input KS Green functions G s . 

Regarding the energy expression of Eq. [7] as a func- 
tional of G yields the well-known Klein functional (24| . 
This functional is stationary (i.e., SE[G)/8G = 0) at the 
self-consistent G of the Dyson equation • It has fur- 
ther been shown 2^, 25 1 that evaluating the Klein func- 
tional (using the GW approximation for $) with the KS 
reference Green function G s one obtains the RPA total 
energy in Eq.[5](see also the supplemental material [2lj|). 
This offers a second way to look at the difference between 
sc-RPA and sc-GW: the sc-RPA energy corresponds to a 
mininum of the Klein functional within a variational sub- 
space of non-interacing KS Green functions, whereas the 
sc-GW total energy corresponds to a stationary point of 
the Klein functional in a larger variational space includ- 
ing both nonintcracting and interacting Green functions. 
However, we emphasize that this stationary point is not 
necessarily a minimum fiol.1241]. In practical calculations, 
the sc-GW total energy is actually above the sc-RPA en- 
ergy as we will show in this Letter. 

For a quantitative comparison between sc-GW and sc- 
RPA, we choose the Galitskii-Migdal (GM) formula ^ 
for the computation of the sc-GW total energy, which 
yields the same result as the Klein functional with the 
self-consistent Green function input, but whose evalua- 
tion is simpler 27J, |28[ . The GM formula can be expressed 
as [HI: 



e gw [ G ] = T + £ cxt + En + E x 



E 



GW 



(10) 



where all terms on the right hand side of Eq. [10] are re- 
garded as functionals of the Green function G. T is the 
kinetic energy of the interacting system and Ef w the 
so-called GW correlation energy defined as 



E, 



GW 



G] 



dw 
2^ 



Tr{GM£ c ^[G]M} 



111) 



Here T,f w is the correlation part of the GW self-energy. 
The evaluation of Eq. [10] with a HF (PBE) Green func- 
tion is referred to as G W @HF (G W @PBE) total en- 
ergy. 

An inspection of Eq. 2] and [11] reveals that the differ- 
ence in the sc-GW correlation energy and the sc-RPA 
correlation energy is twofold. First, the sc-GW expres- 
sion is evaluated with an interacting Green function as 
opposed to a Kohn-Sham one in sc-RPA. Second, the 
kinetic correlation energy - i.e., the difference between 
the full kinetic energy and that of the non-interacting 
KS system - are included in Eq. [4] through the coupling 
constant integration, whereas in sc-GW the correlation 
term is purely Coulombic. To facilitate a term-by-term 
comparison between sc-GW and sc-RPA total energies, 
we separate £? RPA into the Coulomb correlation energy 
C/ RPA and the kinetic correlation energy T RPA : 



U, 



RPA 



— Tr 

2tt 



E? W [G B ] 



(12) 



and T RPA = E* FA -U]F A . The kinetic energy in sc-RPA 
is then given by: 



rpRPA _ y 1 rpRPA 



(13) 



With this reorganization of terms, the kinetic energy in 
sc-GW can be directly compared to T RPA , and similarly 
Ef w to Z7 RPA . Now the only factor responsible for the 
difference in these different pairs of terms arises from the 
difference in the input Green functions used to evaluate 
them. 

We determined the RPA correlation potential follow- 
ing the direct minimization scheme of Yang et al. [29| . 
The resulting orbitals and eigenvalues were used to eval- 
uate the sc-RPA total energy from Eq. [5] We refer to a 
previous publication for details of the sc-RPA implemen- 
tation [30(. The sc-GW method - based on the iterative 
solution of Eqs. [3] and [6] at A = 1 - has been implemented 
in the all-electron localized basis code FHI-aims [3l| , as 
explained in more detail in Refs. 28, [32|. The sc-GW 



total energy was then obtained from Eq. 1101 

We now turn to an assessment of sc-RPA and sc-GW 
for the potential energy curve of H2. In the following, 
we explicitly refer to non-self-consistent calculations by 
appending the suffix @input to label the Green function 
used as input. Figure [2] reports the total energy of H2 for 
different flavors of GW and RPA. For comparison we re- 
produce the full configuration interaction (CI) potential 
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Figure 3. (Color online) Left panel: difference between the sc- 
GW and sc-RPA total energy (AE to t), the correlation energy 
(AE C = E° w - t/ c RPA ) and the remaining terms (A(E to t - 
E c )). Right panel: breakdown of the remaining term into the 
difference of the Hartree (AEu), the external (A-E e xt), the 
exchange (A_E X ), and the kinetic energy (AT). 



energy curve of H2 [18|, that provides an exact refer- 
ence for this system. We also report the total energy of 
H2 evaluated from a beyond-GW/RPA approach that in- 
corporates second-order screened exchange (SOSEX) and 
renormalized single-excitations in the self-energy [33j , re- 
ferred to in the following as renormalized second order 
perturbation theory (rPT2) 34 1. As reported previ- 



ously 1|| 35H371]. non-self-consistent RPA overestimates 
the total energy of H 2 at the equilibrium bond length. 
Around the equilibrium distance, the RPA total energy 
based on exact exchange (OEPx) and sc-RPA arc almost 
identical and overestimate the total energy by approxi- 
mately 0.8 eV, compared to full-CI. At intermediate bond 
distances and in the dissociation region we see a lower- 
ing of the sc-RPA energy compared to RPA@OEPx. The 
spurious "bump" [l6l . |35| , present in all RPA calculations 
for H2 and other covalently bonded molecules, is reduced 
in sc-RPA but is still present. The total energy stays be- 
low the full-CI energy throughout, indicating a general 
overestimation of the bonding and dissociation regions. 



In agreement with Stan et al. [38[, sc-GW provides an 
accurate total energy for H2 close to equilibrium. For the 
Galitskii-Migdal framework, self-consistency is crucial as 
G W @HF and G Wo@PBE largely overestimate the to- 
tal energy. In contrast, the Klein functional evaluated 
with the HF Green function (RPAQHF) provides similar 
results to sc-GW. sc-RPA and sc-GW thus provide a 
qualitatively similar description of the energetics of the 
covalent bond of H2, which results in a slight overestima- 
tion of the total energy (see Table [l|. However, sc-GW is 
in better agreement with full-CI. Most interestingly, the 
sc-GW energy is higher than the sc-RPA one. This is 
in contrast to the exchange only case, in which the HF 
total energy is always lower than (or equal for a two elec- 
tron system) the OEPx energy |l2j ■ This is expected, as 
HF is variational and the local potential in OEPx pro- 
vides an additional constraint that increases the energy. 
Conversely, the total energy in sc-GW has to be higher 



than in sc-RPA, because the variational procedure yields 
a maximum at the self-consistent Green function jlOl |24 1 . 

In the dissociation region, sc-RPA and sc-GW deviate 
markedly. For sc-RPA the dissociation energy is below 
the full-CI energy and is in rather good agreement with 
the reference. sc-GW, on the other hand, fails dramati- 
cally in the dissociation limit and with 24.5 eV underes- 
timates the total energy considerably. On the plus side, 
sc-GW dissociates monotonically and therefore does not 
show the unphysical "bump" present in all RPA-based 
approaches. Again, both non-self-consistent GoW)@HF 
and GoW)@PBE energies give better agreement with the 
reference curve than sc-GW. 

One could surmise that this qualitatively different be- 
havior originates from the different treatment of the ki- 
netic energies that we discussed earlier. Figure O how- 
ever, shows that this is not the case. At equilibrium 
the kinetic energy in sc-GW differs only slightly from 
the sc-RPA kinetic energy defined in Eq. [13j This in- 
dicates that in the bonding regime the AC framework 
correctly reproduces the kinetic energy of an interact- 
ing system. At larger bond distances, the kinetic en- 
ergy differs increasingly in the two approaches. How- 
ever, this effect is averaged out by an opposing change 
in the external energy that arises from an increasing de- 
viation in the electron densities. The same is observed 
for the Hartree and the exchange energy, although the 
absolute magnitude of the effect is smaller. The to- 
tal energy difference between sc-GW and sc-RPA can 
be finally ascribed to the Coulomb correlation energy 
AE C = E? w - t/ c RPA = E^ W [G] - E^ W [G S ], as the left 
panel of Fig. [3] demonstrates. Close to equilibrium AE C 
is of the order of 1 eV, but increases to approximately 4 
eV at larger bond lengths. This illustrates that it mat- 
ters decisively whether the correlation energy is evalu- 
ated with the interacting sc-GW or the non-interacting 
sc-RPA Green function. 

Why the difference is so pronounced at dissociation 
is still an open question. A potential explanation can be 
found in the inverse dependence of the RPA Coulomb cor- 
relation energy on the gap between the highest occupied 
and the lowest unoccupied molecular orbital (HOMO and 
LUMO, respectively). This is exemplified by the right 
panel of Fig. [H which shows the inverse of U c as a func- 
tion of the gap for a simplified two level system. The large 
value of U c obtained from sc-RPA for H2 at dissociation 
can therefore be traced back to the small HOMO-LUMO 
gap (left panel of Fig. 4) of the RPA Green function. In 
contrast, due to the spatial non-locality of the self-energy, 
the HOMO-LUMO gap of the HF and GW Green func- 
tions is much larger at every given bond distance. This 
leads in turn to a smaller Coulomb correlation energy for 
sc-GW and HF-based perturbative methods. 

In conclusion, we have compared MBPT in the GW 
approximation to DFT in the RPA. We found that the 
density functional description is superior at dissociation, 
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Figure 4. (Color online) Left panel: HOMO-LUMO gap ex- 
tracted from the sc-GW spectral function and from HF, PBE 
, and sc-RPA eigenvalues. The HOMO-LUMO gap evalu- 
ated from PBE and HF total energy differences (ASCF) is 
included for comparison. Full configuration interaction (CI) 
calculation were done with a aug-cc-pVDZ basis set. Right 
panel: Inverse of the RPA Coulomb correlation energy [/ C RPA 
for a two-level model as a function of the HOMO-LUMO gap. 



yielding a total energy in qualitative agreement with the 
exact energy along the entire dissociation curve. These 
results illustrate how MBPT and DFT based approaches 
deal with multi-reference ground-states. Kohn-Sham 
DFT is a single determinant approach, which for the ex- 
act exchange-correlation functional would yield the exact 
energy. At dissociation the single particle states in this 
determinant become degenerate as the Kohn-Sham gap 
closes, which mimics the multi-dctcrminantal solution 
found in full CI. The same effect in Green function the- 
ory has to be achieved by the right (potentially infinite) 
set of diagrams. This has profound consequences for the 
development of new computational electronic structure 
schemes and highlights the importance of going beyond 
GW (vertex corrections) to account for those effects. Our 
results for bond breaking/formation of H2 illustrate that 
RPA-based DFT is a strong contender for a universally 
applicable electronic-structure theory. 
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